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Abstract 

Electromagnetic effects are increasingly being accounted for in lattice quantum chro- 
modynamics computations. Because of their long-range nature, they lead to large finite- 
size effects over which it is important to gain analytical control. Nonrelativistic effective 
field theories provide an efficient tool to describe these effects. Here we argue that some 
care has to be taken when applying these methods to quantum electrodynamics in a finite 
volume. 


State-of-the-art lattice quantum chromodynamics (QCD) computations have reached such 
a high level of accuracy (see e.g. |1 and references therein) that small electromagnetic correc¬ 
tions, and other isospin breaking effects, are becoming important. These effects will have to be 
accounted for more and more systematically if the results of lattice calculations are to continue 
to be used to test the standard model and to search for new physics in increasingly precise 
experiments. Moreover, electromagnetic effects in hadrons are important in themselves. For 
instance, they are critical for understanding Big Bang nucleosynthesis as well as many proper¬ 
ties of atomic nuclei or for determining the up and down quark masses. As a result, increasing 
attention is being focussed on including quantum electrodynamics (QED) corrections in lat¬ 
tice QCD calculations. A number of results concerning hadron and quark masses have been 
obtained, in the electroquenched approximation 2j-‘8) and in full QCD+QED [9j-[l2 . In ad¬ 
dition, a method for including QED corrections in the lattice calculation of hadronic matrix 
elements has been proposed in TTl) . 

A very important issue in such calculations is finite-volume effects.^ Indeed, the vanishing 
mass of the photon implies that the leading finite-volume corrections are proportional to 

1 Throughout this paper we will be concerned with situations in which the linear dimensions of the lattice 

are much larger than the Compton wavelengths and the internal-structure length-scales of the particles under 

consideration. 
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inverse powers of the spatial dimension, L, of the lattice, instead of being exponentially 
suppressed, as they are in pure QCD calculations. These power corrections represent a large 
fraction-tens of percent-of the computed electromagnetic effects for typical lattice sizes L ~ 
3fm, as shown in j5j[8,12 , where controlled infinite-volume extrapolations were performed. 
These corrections must be appropriately subtracted to obtain accurate results. Thus, it is 
important to know their precise analytical form and to constrain them as much as possible. [^] 
An elegant and efficient method to determine the functional form of these effects was pro¬ 
posed and worked out in 14 . It is based on nonrelativistic effective field theories (NREFTs) 


and allows to compute these corrections in a systematic expansion in powers of 1/L for any 
spin-0 or 1/2 massive particles with (or without) internal structure, such as hadrons or nuclei. 
In this approach, the corrections are determined by the mass and by well-defined electromag¬ 
netic properties of the particle, such as its charge, charge radius, magnetic moment, etc. 
However, as pointed out in 12 and explained in more detail below, the results obtained 
in 14 do not fully agree with those of point-particle calculations performed in 12 . Since 
many checks of the point-particle calculation were performed in [12] and both results cannot 


be correct, we have been led to investigate the possible source of the discrepancy. As explained 
below, this investigation has led us to uncover a subtlety in the application of NREFTs to 
finite-volume QED, which goes beyond the calculation discussed in the present paper. [^] 

To expose the problem and its solution, we focus on the calculation of the finite-volume 
corrections to the pole mass of spin-1/2 particles at O(a). 


As 


m 


12, 14 


we work with 

periodic boundary conditions in the QED^ formulation of QED in a finite volume. In this 
formulation, momentum modes of the photon field with k = 0, for all values of k$, are 
eliminated from the theory. This approach was first proposed in 1161 and, as shown in 12 


has many theoretical and practical advantages. Here we work with an infinite time direction, 
as it simplifies analytical calculations and is equivalent to the finite time extent, T, formulation 
with periodic boundary conditions, up to corrections which are smaller than any inverse power 
in T 
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Finite-volume, electromagnetic corrections to the mass of a particle can be obtained from 
the difference of its on-shell, electromagnetic self-energy at rest, in finite and infinite volumes. 
Physically, the corrections in inverse powers of L arise mainly from the particle under consider¬ 
ation exchanging a photon with itself around the periodic three-volume. In |14j, Davoudi and 
Savage (DS) compute this effect for a generic spin-1/2 particle using NREFT to N 3 LO, i.e. to 
order 1/L 4 . The use of NREFT is entirely justified here, because we are after an asymptotic 
expansion of particle properties in inverse powers of L. NREFT is a low-momentum effective 
theory of QED and in QED L , it provides expansions of low-momentum particle properties in 
powers of the infrared momentum cutoff, 2i r/L. 

To understand the NREFT approach, we have repeated their calculation. We stop at 
N 2 LO, i.e. up to and including terms proportional to 1/L 3 , because this is the order at which 
the subtlety, described below, first arises. We use the same NREFT lagrangian as they do, 


2 The two leading finite-volume corrections, proportional to 1/L and 1/L 2 , are given in terms only of 


the particle’s charge and infinite-volume mass [12 14 . This important feature allows precise infinite-volume 


extrapolations of the QED contributions to particle masses in lattice calculations. In particular, the finite- 
volume corrections do not depend on spin nor on particle structure |12j[f4]. As shown in [12], this universality 
follows, under very general hypotheses, from QED Ward-Takahashi identities and the work of Liischer on the 
analyticity properties of propagators and vertex functions 


15 


3 It is important to note that this subtlety does not affect the two leading, universal, finite-volume corrections 
mentioned above, as it only enters at order 1/L 3 , as shown below. 
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Figure 1: NREFT self-energy diagrams which contribute to finite-volume effects on the mass 
of a spin-1/2 particle, up to order 1/L 3 . Their contributions were taken into account in fifl. 


The solid line corresponds to the fermion propagator; the dashed line to the temporal photon 
propagator; the wavy line to the transverse photon propagator; the cross to a kinetic term 
insertion; the simple vertices with these lines, to the fermion-temporal-photon interaction; to 
the fermion-transverse-photon interaction; to the fermion-to-two-transverse-photons interac¬ 
tion; the square vertex to the Fermi interaction; the crossed vertex to the Darwin interaction; 
the dot vertex to the spin-orbit interaction. Under each diagram we indicate its contribution 
to the finite-volume correction on the fermion mass. Calcidations are performed in Coulomb 
gauge. An explicit factor of 2 appears in the single Fermi, Darwin and spin-orbit vertex terms. 
This indicates that to each of the three diagrams shown corresponds one in which the vertices 
are switched. These new diagrams give the same contributions as the original diagrams. The 
sum of all the contributions shown in the figure yields the result in Eq. 0- 
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where ^ is a two-component spinor which annihilates a particle, D^ + ieQA fl , cf = 

Q + k + O(a), cd = Q + 4m 2 (r 2 )/3 + O(a), eQ is the charge of the particle, m its infinite- 
volume mass, (r 2 ) its mean-squared charge radius and k its anomalous magnetic moment. 
The diagrams which contribute up to N 2 LO are shown in Fig. [lj They yield the following 
hnite-volume corrections: 


Ci 
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is 


where m(L) is the value of the particle’s mass in QEDp and Ci = —2.837297(1) 12,25 -27 
a known constant. This result fully agrees with DS’s. 

However, as noted above and in jl2|, if we reduce the corrections of Eq. © to the point 
particle case by setting (r 2 ) = 0 and n = 0, we find that it disagrees with the one that we 
computed directly in spinor QED. Indeed, in the reduced corrections, the term proportional 
to (irctQ 2 /m 2 L 3 ) has a coefficient of 3/2 instead of 3, as found in 


12 


Moreover, in 


12 


we 

performed a precise, dedicated numerical study of finite-volume effects on the pole mass of 
a point-like fermion in QEDp. The study was conducted with fixed bare parameters and six 
lattice sizes ranging from L/a = 24 to L/a = 128, with a the lattice spacing. A fit of the 
measured, finite-volume mass, to a polynomial in a/L, highly favored our value for the 1/L 3 
coefficient over the one in |14j. To explain where this discrepancy comes from, we must now 
describe the main features of our relativistic calculation 1121. 

The computation in spinor QED is a straightforward asymptotic expansion in large L, 
of the difference between the finite and infinite-volume, on-shell self-energies of the parti¬ 
cle 12 . As mentioned above, we consider this difference at O(a), i.e. we consider the usual, 
one-loop sunset diagram in which a photon is emitted and re-absorbed by the particle. The 
asymptotic expansion is most straightforwardly performed using a Poisson summation for¬ 
mula. In that approach, the corrections in powers of 1/L result from the nonanalyticities 
in the integrand/summand, associated with intermediate states going on shell in the domain 
of integration. The obvious singularities that arise in the present case are the particle and 
the positive and negative energy photon poles. The antiparticle pole, which is 2m away in 
energy, is only expected to contribute terms which fall off exponentially with 2mL. This is 
an illustration of the decoupling of antiparticle modes that leads to the NREFT for a single, 
massive particle, used in |l4 . 

Upon closer inspection, however, we find that this expectation is incorrect. Analyzing the 
different contributions shows that the antiparticle pole contributes to the contentious term 
of order l/(m 2 T 3 ), even though its propagator cannot go on shell for the given kinematics. 
This surprising result is due to the fact that the contribution of the zero modes of the photon 
are omitted from the loop sum. At order 1/L 3 , the finite-volume corrections come not only 
from singularities of the summand/integrand in the domain of integration, but also from 
the explicit subtraction of the photon zero modes. Since these modes couple particles to 
antiparticles, the latter also play a role in the calculation of finite-volume effects. 
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Figure 2: Antiparticle contribution to the self-energy of a spin-1/2 particle which arises from, 
the four-fermion lagrangian of Eq. 0- The double line corresponds to the antiparticle prop¬ 
agator. The vertex corresponds to the vector four-fermion coupling. Under the diagram we 
indicate its contribution to the finite-volume corrections on the fermion mass. The diagram 
contributes at order 1/L 3 and provides the term which is missing in the calculation of [If]. 


In the language of NREFT, we reach the same conclusions if we explicitly include an¬ 
tiparticle degrees of freedom, in addition to the usual contributions of antiparticles to higher- 
dimension particle operators. Thus, to the lagrangian density, for particle fields given in 
Eq. @ and used in |14], we add the one for antiparticle fields, Xi and include direct couplings 
of particles to antiparticles. Then, the lagrangian density becomes: 

£n 2 lo = A/ 1 + Cx + £ 4 / + 0(j l 4 ) , (3) 

where x is a two-component spinor which annihilates antiparticles. In Eq. ([ 3 ]), C x is obtained 
from C with the replacements if — > x and Q —> —Q. The relevant four-fermion interaction, 
at 0(a), originates from particle-antiparticle annihilation in the triplet channel and is given 

by [28] 

C Af = d v -°^(i[ ji a(j 2 X*) ■ (x T V 2 v/>) + 0(a 2 ,p 4 ) , (4) 

where the cr,, i = 1,2, 3, are the Pauli matrices. For point particles, dy = —TrQ 2 +0(a) 0® 
From this four-fermion lagrangian, it is straightforward to compute the contribution of 
the antiparticle to the finite-volume effects in the particle’s mass. It is given by the self-energy 
diagram of Fig. [2] We find this contribution to be 

A 4 ML) = Y, 1 > ( 5 ) 

frO 


where q is the momentum of the antiparticle in the loop and the sum, represents the 

difference between the sum over the finite-volume modes and the infinite-volume integral, i.e. 


1 \] 1 \ ^ _ f d 3 q 
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fro 
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We now argue that the q = 0 modes must be eliminated from the finite-volume sum, even 
though it is, initially, only photon zero modes which are removed. This is where the subtlety 
enters in NREFT. If the q = 0 antiparticle modes were present, as one might guess, then the 
finite-volume correction of Eq. ([ 5 ]) would vanish and antiparticle degrees of freedom would 
not contribute. However, one must remember where internal particle or antiparticle lines 
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in NREFT come from. In relativistic QED, an internal particle or antiparticle line in a 
diagram, such as the self-energy diagram under consideration, is produced at a vertex with a 
photon. But in QED^ that photon cannot have a vanishing three-momentum. Therefore in 
such diagrams, where all external particle lines have vanishing three-momenta, the internal 
antiparticle lines cannot have vanishing three-momenta. This justifies the omission of the 
q = 0 antiparticle modes in the contribution of Eq. ([H]). 

Now o 1 = — 1) so the full NREFT expression for the finite-volume correction to the 
mass m of a spin-1/2 particle of charge Q is, to 0(1/L 3 ) in the presence of electromagnetism, 
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Using the value of dy for a point particle given after Eq. Q, we find that this additional 
correction adds to the ones in A DS m(L) exactly the ( 3iraQ T /2m 2 L 3 ) term which is missing 
to reproduce the relativistic point-particle result of |12|. 

It is worth noting that the result of Eq. 0 can also be obtained by only reinstating 
the q = 0 modes of the antiparticle field in the finite-volume NREFT. Then one uses the 
corresponding terms in the lagrangian to subtract these modes’ contribution from the finite- 
volume self-energy. In that approach, the 1 that appears in Eq. (5) would directly be 

replaced by —1. One need consider neither the spatial modes of the antiparticle field in finite 
volume, nor any antiparticle modes in infinite volume. 

To conclude, we have considered the calculation of finite-volume corrections to the pole 
mass of a charged spin-1/2 particle, in the presence of electromagnetism. We have explained 
how the NREFT calculation of [14] can be reconciled with the relativistic QED result of 12 


In the process, we have shown that there are subtleties associated with applying NREFTs to 
the calculation of finite-volume effects in QED. In particular, we have argued that antiparti¬ 
cle degrees of freedom must be dealt with carefully, because photon zero modes are treated 
differently in finite and infinite volumes. Indeed, those modes can couple antiparticles to 
particles. Therefore, NREFT calculations of particle properties must account for the contri¬ 
bution of antiparticles to the removal of the photon zero modes in finite volume. However, 
once this contribution is suitably accounted for, the NREFT approach constitutes an elegant 
and efficient formalism to calculate finite-volume corrections. 
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